A SINGULAR LIMIT FOR COMPRESSIBLE 
ROTATING FLUIDS 



EDUARD FEIREISL, ISABELLE GALLAGHER, AND ANTONIN NOVOTNY 



Abstract. We consider a singular limit problem for the Navier- 
Stokes system of a rotating compressible fluid, where the Rossby 
and Mach numbers tend simultaneously to zero. The limit problem 
is identified as the 2-D Navier-Stokes system in the "horizontal" 
variables containing an extra term that accounts for compressibility 
in the original system. 

1. Introduction 
Consider a scaled Navier-Stokes system in the form 

(1.1) d t g + div x (gu) = 0, 

(1.2) d t (gu) + div^u g> u) + ~(g x gu) + \v x p(g) = div x S(V x u), 
with the viscous stress tensor 

(1.3) S(V s u) ^//(v^u + V^u-^div^ul), /i>0, 
and 

g= [0,0,1]. 

Here g = g(t,x) > denotes the density and u(t,x) = [wi, w 2 , «3](t, a;) 
denotes the velocity of the fluid. Problem (11.11 - II. 2p arises in meteo- 
rological applications, modeling rotating compressible fluids with the 
rotation axis determined by g and the Rossby and Mach number pro- 
portional to a small parameter e. 

We consider a very simple geometry of the underlying physical space, 
namely an infinite slab Q bounded above and below by two parallel 
planes, 

(1.4) fi = M 2 x(0,l). 



The work of E.F. was supported by Grant 201/08/0315 of GA CR as a part of 
the general research programme of the Academy of Sciences of the Czech Republic, 
Institutional Research Plan AV0Z10190503. 

The work of I.G is partially supported by the French Ministry of Research grant 
ANR-08-BLAN-0301-01. 

1 



2 EDUARD FEIREISL, ISABELLE GALLAGHER, AND ANTONIN NOVOTNY 

The velocity u satisfies the complete slip boundary conditions, 

(1.5) u-n = M 3 |an = 0, [Sn] x n\ d n = [£2,3, -#1,3, 0] | 9 n = 0. 
For the initial data 

(1.6) q(0, •) = g 0)£ , u(0, •) = Uq )6 , 

our goal is to study the asymptotic behavior of the corresponding solu- 
tions g e , u £ for e — y 0. We focus on the interplay between the Coriolis 
force, here proportional to a singular parameter 1/e, and the acous- 
tic waves created in the low Mach number regime. In particular, we 
neglect: 

• stratification due to the presence of gravitation, here assumed 
in equilibrium with the centrifugal force; accordingly the action 
of the centrifugal force is also neglected; 

• the effect of a boundary layer (Ekman layer), here eliminated 
by the choice of the complete slip boundary conditions. 

We consider ill-prepared initial data, specifically, 

( Qo,e = f? + er ,e, with {r , £ } £ >o bounded in L 2 D L°°(Q), 



;i-7) 



for some positive constant g, 



{u , £ } £ > bounded in L 2 n L°°(fi; R 3 ). 

Because of the prominent role of the "vertical" direction g in the 
problem, we introduce the "horizontal" component = [fi,f2,0] of 
a vector field v, together with the corresponding differential opera- 
tors Vfc, div h , and, notably, curl ft , which is represented by the scalar 
field 

curl/Jv] = d Xl v 2 - d X2 v x . 

Let g E , u £ be a solution of problem fll.ltll.6l) . Introducing a new quan- 
tity 

Qe ~ Q 



which satisfies 



d t r £ + -i-diy^u) + div x (r £ u) = 0, 



we easily check that if 

r £ — > r, u £ — > U in some sense, 
then, at least formally, the limits satisfy a diagnostic equation 

(1.8) gxU + ^V/^, 

Q 
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which in turn implies that 

(1.9) r = r(zi, x 2 ), U = [V h , 0], U fc = U h (xi, x 2 ). 

Moreover, as we will see below, div^U = div^U/i = 0, and denot- 
ing V^r the vector (d X2 r, —d Xl r), 

(1.10) d t (A h r - -^-r) + Vir ■ V h {A h r) = ^A\r. 

Thus r may be interpreted as a stream function associated to the vec- 
tor field XJh, therefore fll.lOp can be viewed as a 2D Navier-Stokes sys- 
tem describing the motion of an incompressible fluid in the horizontal 
plane M 2 , supplemented with an extra term (1 / 'p' (g))d t r . 

The main goal of the present paper is to provide a rigorous justifi- 
cation of the target system fll.lOp in the framework of weak solutions 
to the primitive equations (II. ip . (II. 2ft . In Section EJ we introduce the 
weak solutions to both systems, recall their basic properties, and state 
our main result. In Section [3j we derive the necessary uniform bounds 
on the family of solutions {o e ,u e } e>0 , and pass formally to the limit 
when e — > 0. In Section HI the associated wave equation describing 
propagation of the acoustic waves in the low Mach number regime is 
introduced. Using the celebrated RAGE theorem, we show that the 
acoustic energy tends to zero, at least locally in space. The proof of 
the main result is completed in Section |5j 

2. Preliminaries 

To begin, we point out that system (11.11 - II. 3p . endowed with the 
boundary conditions (II. 5p can be recast as a purely periodic problem 
with respect to the vertical coordinate x 3 provided g, u\, u 2 were ex- 
tended as even functions in the x%— variable defined on 

while W3 is extended to be odd in x% on the same set. A similar con- 
vention is adopted for the initial data. 

2.1. Weak solutions. We shall say that functions g, u represent a 
weak solution to problem (11.11 - [L6|) in (0,T) x Q if: 

• g > 0, g e L°°{0, T- L^{Q)) for a certain 7 > 3/2, 

u6L 2 (0,T;r> 2 (n;l 3 )); 

• equation of continuity (11.11) is satisfied in the sense of renormal- 
ized solutions, namely 

J J ( K (g+K0))dt^+(g+K0)) u -^xV+(K0)-b\g)g)div x u^ dx dt 
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(eo,e + 6(eo, e )V(o, •) dx 

for any b G C°°[0, oo), b' G C£°[0,oo), and any test function 
<peC?([0,T) xfi); 
• p — p(g) G ^((O, T)xO), momentum equation (II .2p is replaced 
by a family of integral identities 



(2.2) 




gu- d t f + q(u®u) : V x ^ + -(gx gu)-<f + —p(g)div x <fj dx dt 




(2.3) 



/ §(V x u) : V x p dx dt - / £o, e u , e • (p(0, •) dx 
o Jn 

for any ^ G C 6 °°([0,T) x ft; 
» the energy inequality 

J (^g\u\ 2 + ^E(g,g)^j (r,-) dx + J S(V x u) : V x u dx dt 



Q^0, £ |u , £ | 2 + ^^(^O.e, ^)^) dx 

holds for a.a. r G (0, T), where 

E(q,q) = H(g) - H'{q){q -q)- H(g), 

with 

H( e ) = q J' dz. 

Note that, by virtue of hypothesis (|1.7p . the quantity on the right- 
hand side of (12. 3 p is bounded uniformly for e — > 0. 

Existence of global-in-time weak solutions to problem (11.11 - II. 6p 
can be established by the method developed by P.-L. Lions (TD], with 
the necessary modifications introduced in [S] in order to accommodate 
a larger class of physically relevant pressure-density state equations, 
specifically, 
(2.4) 

p G C^O, oo), p(0) = 0, pYg) > for e > 0, lim ^4 = Poo > 0, 




for a certain 7 > 3/2. 
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2.2. Main result. The main result of the present paper can be stated 
as follows. 

Theorem 2.1. Assume that the pressure p satisfies \2.J$ . 

Let {g E , u £ } e>0 be a family of weak solutions to problem U.l\ - \1.6\) 
in (0, T) x Q, where Q is specified through with the initial data 

satisfying ( [1. 7| ) ; where 

r 0j£ — > r weakly in L 2 (fi), u 0ie — > Uo weakly in L 2 (f2;M 3 ). 

Then after taking a subsequence, the following results hold 

r e = r weakly-(*) in L°°(0, T; L 2 {Q) + L 7 (ft)), 

u e ->• U weakly in L 2 (0, T; W 1,2 (fi; K 3 )), 

and u £ ->■ U strongly in L 2 oc ((0,T) x fi;]R 3 ), 

where r and U satisfy U.8\) . div x U = ; and, moreover, the stream 
function r solves equation /11.10\) in the sense of distributions, supple- 
mented with the initial datum 

(2.5) r(0,-) = r, 

where r e Fy 1,2 (R 2 ) zs t/ie unique solution of 

-A h r + — ]—r = g [ curUU ,h dx 3 + / r dx 3 . 
P{q) Jo Jo 

If, in addition, curlhUo^ € L 2 (Q), then the solution r of U.10\) is 
uniquely determined by Ii2.5\) and the convergence holds for the whole 
sequence of solutions. 

The remaining part of the paper is devoted to the proof of The- 
orem 12.11 The crucial point of the proof is, of course, the strong 
(a. a. pointwise) convergence of the velocity field that enables us to 
carry out the limit in the convective term. Here, the desired point- 
wise convergence will follow from the celebrated RAGE theorem, to- 
gether with the fact that the wave propagator in the associated acoustic 
equation commutes with the Fourier transform in both the horizontal 
variables (x\,X2) and the vertical variable 23. 

2.3. Related results. This work is a contribution to a general re- 
search direction consisting in studying singular limits in PDEs aris- 
ing in fluid mechanics. Without giving an extensive bibliography, one 
should refer for the first works in this line to Klainerman and Majda 
[9] and Ukai [11] for the incompressible limit (actually [11] is probably 
the first work in which dispersive estimates were established in order to 
prove strong convergence in the whole space), followed by Desjardins et 
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al. [5] and [B]. In the context of rotating fluids one should mention the 
important work of Babin, Mahalov and Nicolaenko [TJ, as well as the 
book [3] and references therein; one also refers to [7\ for a survey. Few 
studies combine both rotation and compressible effects. We refer to 
Bresch, Desjardins and Gerard- Varet [2] for an analysis in a cylinder, 
where the well prepared case is studied precisely; the ill prepared case 
is also addressed but only a conditional result is proved. 

3. Uniform bounds 

We start reviewing rather standard uniform bounds that follow di- 
rectly from the energy inequality (12. 3p . To this end, it is convenient to 
introduce a decomposition 

h = [h] css + [h} rcs , where [h] css = ip(g e )h, 

ip G C^°(0, oo), < tj) < 1, tj) = 1 in a neighborhood of ~g 
for any function h defined on (0, T) x Q. It is understood that the es- 
sential part [ft.] ess is the crucial quantity that determines the asymptotic 
behavior of the system while the residual component [h] res "disappears" 
in the limit e — y 0. 

As already pointed out, our choice of the initial data ( I1.7P guaran- 
tees that the right-hand side of energy inequality (12.31) remains bounded 
for e — > 0. After a straightforward manipulation, we deduce the fol- 
lowing estimates: 

(3.1) {^Q~sU £ } £>0 bounded in L°°(0, T; L 2 (fi; M 3 )), 

(3.2) {[r £ ]ess} £ >o bounded in L°°(0, T; L 2 (fi)), 

(3.3) ess sup ||[^]res|ll 7( m < £ 2 c, 

te(o,T) 

(3.4) ess sup ||[l] res ||n(n) < £ 2 , 

te(o,T) 

and 
(3.5) 

i V x u £ + V^u £ - -div x u £ ll bounded in L 2 ((0,T) x Q;R 3x3 ). 

I 3 J £>0 

In addition, it is easy to observe that (13. 2p . (13. 3p yield 

(3.6) q £ -+-Qm L°°(0,T;L 7 + L 2 (fi)), 

which, together with (13. ip . (13 .6p and the standard Korn inequality, 
gives rise to 

(3.7) {u e } £>0 bounded in L 2 (0, T; W 1 ' 2 ^; M 3 )). 
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In accordance with (13. 2D . (13 .3D . we may assume that 

(3.8) [r £ ] css -> r weakly-(*) in L°°(0, T; L 2 (Q)), 
and, taking H3.4D into account, 

(3.9) [r £ ] res -> in L°°(0, T; L 9 (fi)) for any 1 < q < min{ 7 , 2}. 
Moreover, by virtue of (13.71) . 

(3.10) u e -»■ U weakly in L 2 (0, T; W 1,2 (Q; M 3 )), 

passing to suitable subsequences as the case may be. 

Letting e — > in the weak formulation of the continuity equation 
(12.11) . with b = 0, we obtain 

(3.11) div x U = a.a. in (0,T) x Q. 

Finally, multiplying momentum balance (12. 2p by e, we recover (II. 8p 

-E/ 2 

in particular, r = r(x\ 1 X2) is independent of the vertical variable, and 

(3.13) U h = U h (xi,x 2 ), div h U h = 0, 

which, together with (13. lip , implies U3 is independent of X3. However, 
as U satisfies the complete-slip boundary conditions (II. 5p on dfl, we 
may infer that 

(3.14) U 3 = 0. 

4. Propagation of acoustic waves 

Assume from now on, to simplify notation, that p'(j)) = 1. Sys- 
tem (12. ip . (12. 2p can be written in the form 

(4.1) ed t r £ + div^V, = 0, 

(4.2) ed t V £ + (g x V £ + V^) = 
where we have set 



(3.12) 



p'(g)V h r, d 3 r = 0, 



Qe-Q , r 

r £ = , V e = g £ u £ , 



and 



f e = diY x S(V x u E )-diY x (g e u E <S>u E )-—V x ^p(g e )-p'(g)(g E -g)-p(g) 

More precisely, system (14. ip . (14. 2 p should be understood in the weak 
sense: 

(4.3) / I (er £ d t <p + V £ • V x ip) dx dt = -e [ r 0j£ ^(0, •) dx 
Jo Jn v J Jn 
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for any <p E C c °°([0,T) x fi), 
(4.4) 



J J (sV £ -d t (p-(gx V e )-(p+r e dxv x (pj dx dt = -a J <i £ ,ip>dt 



Jn 

for any test function ip G C£°([0, T) x fi; R 3 ), <p ■ n|an = 0, where 
- < i E ,ip >= J (§(V x u £ ) : - {q £ u e <g> u e ) : V^y? 

--^(p(Qe) -p'(q)(Qs -q) - p{gj)dxv x <p) dx. 
It follows from the uniform bounds established in ( 13.11 - 13 . Tj) that 

(4.5) < f £ , p >= y (V] : V x ip + F 2 : V.^) dx, 
with 

(4.6) {^} £>0 bounded in L°°(0, T; L 1 ^; M 3x3 )), 

(4.7) {F 2 } £ > bounded in L 2 (0, T; L 2 (fi; M 3x3 )). 

4.1. Point spectrum of the acoustic propagator. Consider an op- 
erator B defined, formally, in L 2 (f2) x L 2 {Vl; M 3 ), 



B 



r 




V 





div x V 
x V + V x r 



As a matter of fact, it is more convenient to work in the frequency 
space, meaning, we associate to a function v its Fourier transform v 

v = v(£ h ,k), 6, = (£i,6)eK 2 , kez, 



where 



/ / exp ( •x fe ))v(x /l ,X3) dx h exp(-ifcx 3 ) dx 3 

Jo 7r2 v ' 



We investigate the point spectrum of B, meaning, we look for solu- 
tions of the eigenvalue problem 



(4.8) 



div x V = Ar, V x r + g x V = AV, 



or, in the Fourier variables, 



z 

i( J2 + kV *) - Af = °' ^ - -Vi, 0] - AV = 0. 



A SINGULAR LIMIT FOR COMPRESSIBLE ROTATING FLUIDS 



After a bit tedious but straightforward manipulation, we obtain 



(4.9) 



A 



-/i, ii 



1 + K| 2 + k 2 ± v/(l + |£| 2 + P) 2 -4P 



whence the only eigenvalue is A = 0, for which k = 0, and consequently, 
the space of eigenvectors coincides with the null-space of £>, 



(4.10) 



Ker(£) 



[r, V] r = r(xi,X2), 



V = {V l {x 1 ,x 2 ),V 2 {x 1 ,x 2 ),V 3 {x 1 ,x 2 )}, div h V h = 0, V h r={V 2j -V 1 ]Y 

4.2. RAGE theorem. Our goal is to show that the component of the 
field [r e , V e ], orthogonal to the null space Ker(£>) decays to zero on 
any compact subset of Q. To this end, we use the celebrated RAGE 
theorem in the following form (see Cycon et al. [U Theorem 5.8]): 

Theorem 4.1. Let H be a Hilbert space, A : V(A) C H — >■ H a self- 
adjoint operator, C : H — )■ H a compact operator, and P c the orthogonal 
projection onto H c , where 



H = H c ® cl# jspan{u> e H \ w an eigenvector of A}j. 



Then 



exp(-rL4)CP c exp(iL4) dt 



— > for t — > oo. 



C(H) 



In addition to the hypotheses of Theorem I4.1[ suppose that C is 
non-negative and self-adjoint in H. Thus we may write 

1 

T 



exp -i-A Cexp i-A P C X,Y ) dt < h(e)\\X\\ H \\Y\\ H , 



H 



where h(e) — > as e — > 0. Taking Y = P C X we deduce 

r-T / t \ 2 

v^exp ii-A J P C X 
Similarly, for X G L 2 (0,T; H), we have 



(4.11) 



1 

T 



d* <Mc)ll^llSr- 



(4.12) 



1 
1 

T 



t - s 



CP C / exp i A X(s) ds 



T r T 



VCexp ni— i^J P c X(s) 

< h(e) [ exp(-i-A) X(s) 2 ds = h{e) f \\X(s)\\ 2 ds. 
Jo \ e / h J 



L 2 (0,T;H) 
2 

dt ds 

H 

T 
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4.3. Application of RAGE theorem. For a fixed M > 0, we intro- 
duce a Hilbert space 

H = H M = {[r,V] | r(£ h ,k) = 0,V(&,fc) = whenever > M}. 

Let 

P M : L 2 (fi) x L 2 (n;R 3 ) i? M 

denote the associated orthogonal projection onto if a/. 
Our goal is to apply RAGE theorem to the operators 

A = iB, C[v] = P M \xv], x e c^n), o < x < 1, 

considered on the Hilbert space ii^r . 

Going back to system (14.31) . (14.41) . we obtain that 



(4.13) 
where 

and 



dt 



r £ ,M 

V £ ,M 



+ B 



r £ ,M 







[r e ,M, V £ , M ] = -Pjif[r e , V £ 





fe.Af 



G if 



H 



M; 








s 






5 


w 


) -/ 

/ H M JO 



whenever (s,w) G H M . Since 

llwHiy^.^nH/^-^niK 3 ) < C(77l) || w|| W /m+2,2( Q;K 3) < CiVi m+2 || W || L 2( n . R 3), 

we may use the uniform bounds (j4.6j) . (14. 7p in order to conclude that 





< c(M) 



L 2 (0,T;^ M ) 



uniformly for £ — )■ 0. 

Writing solutions to ( I4.13P by means of Duhamel's formula we get 
(4.14) 



V £)M 



exp(iA- 



r £ ,A/(0) 
V £ , M (0) 



+ / exp l A 

'o V £ 





fe,M 



ds; 



whence a direct application of ( 14. lip . ( I4.12p . recalling that the only 
point spectrum is reduced to 0, yields 



(4.15) 



Q 



r £ ,M 

V £)M 



^ in L 2 ((0, T) x K; M 4 )) as e 0, 



for any compact i\" C and any fixed M, where we have denoted 
Q : L 2 (VL) x L 2 (VL-M?) -» Ker(B) 
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the orthogonal projection onto the null space of B. Indeed observe that 

2 



VCQ- 



r £ ,M 

Ve,M 



CQ~ 



r £ ,M 



r £ ,M 
V £i m 



X 



Or 



r £ ,M 
V £ , M 



dx, 



where we have used the fact that Pm and Q commute. 
Finally, a direct inspection of (I4.14p yields 



(4.16) Q 



r £ ,M 



in L 2 ((0,T) x K;R 4 )) as e 0, 



where r and U are the asymptotic limits identified through (13.81 - 1312]) . 



4.4. Strong convergence of the velocity fields. Relations (14. 151) . 
(14. 161) . together with (13.71 - 13. lOf) . may be used to obtain the desired 
conclusion 

(4.17) u e ->■ U in L 2 ((0, T) x K; M 3 ) for any compact K C SI. 
Indeed, by virtue of fl3~9l) . (05jl . fl4~T6|) . we obtain 

P M [u £ ] -> P M [U] in L 2 ((0,T) x if; ]R 3 ) 



for any fixed M, which, together with (13 . 1 0[) and compactness of the 
embedding W l ' 2 {K) ^ L 2 (K), yields f07|) . 

5. The limit system 

5.1. Identifying the limit system. With the convergence estab- 
lished in ( 13.81 - 13.101) , and ( 14. 17[) , it is not difficult to pass to the limit 
in the weak formulation (12.11) . (12.21) . To this end, we take 

y>=[V^,0], ^GL7 c °°([0,T)xfi) 
as a test function in momentum equation (12.21) to obtain 

(5.1) / / (g £ u £ • d t ip + g £ u £ <g> u e : V x ip - -g £ [u £ } h ■ V x ip) dx dt 
Jo Jn v £ J 

= - Qo, £ u 0t£ ■ <p(Q, •) dx + / / S(V x u £ ) : V z y? dx dt. 
Jn Jo Jn 

Moreover, (I4.3|) yields 

(5.2) I [ (r £ dtij + -g £ [u £ ] h ■ V x ^) dx = - [ r ^(0, •) dx. 
Jo Jn v e 7 7n 
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Letting e — > in (I5.ip . (15. 2 p we may infer that 

J J (gUh ■ dtViiP + Q{U h ® U h ) : V^Vjty) + r<9^) dx 

= -J (gV , h • V^(0, •) + r ^(0, ■)) dx 

+ / / nV h U h : V(V^) da; dt. 
Jo Jn 

Moreover, as the limit functions are independent of X3, we get, 
(5.3) / f (g\J h ■ d t Vii> + g\U h ® U ft ] : V h (V^) + ^) dx, dt 

= ~ {I v ° ,h dx3 ) ' v ^ (o ' -) + (1 r ° dx3 ) ^ (o ' '0 d ^ 

+ [ [ tiV h U h : V h (V^) dx h dt 
Jo Jk 2 

for all ip G C£°([0,T) x R 2 ). 

Finally, by virtue of ( 11.81) . = V^r, and (15. 3 j) coincides with a 

weak formulation of (ll.lOp . (I2.5p . We have completed the proof of the 

convergence result, up to a subsequence, of Theorem 12.11 

5.2. Uniqueness for the limit system. In this final section we shall 
prove that the limit system has a unique solution provided the initial 
data are more regular. In order to do so we shall simply write an 
energy- type estimate on the difference of two solutions, called r\ and r 2 , 
associated with two initial data f\ and r<i- This will provide a stability 
estimate, whose immediate consequence will be a uniqueness result. 
Notice that the diagnostic equation (11.81) implies that f should be taken 
in W^iR 2 ). 

The limit system writes 

d t {A h r - r) + Vir ■ V h {A h r) = ^A 2 r 

g 

recalling that for simplicity we have chosen p'ij)) = 1. Multiplying 
(formally) this equation by A^r and integrating over R 2 yields 

|(l|A,r||l 2 + ||V,r||i 2 ) + ^||V,A,r||| 2 =0, 

whence the estimate 

||A,r(t)||| 2 + ||V,r(t)||i 2 + ^ f ||V,A,r(t')||| 2 dt' 

Q Jo 

= \\A h r\\ 2 L2 + ||V fc f ||ia. 
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Now suppose r\ and r 2 are two solutions as described above, and de- 
fine 5 := r 1 — r 2 . Then of course S satisfies 

d t (A h 5 -5) + ■ V h (A h r 2 ) + V^n ■ V h (A h 5) = ^A\8 

Q 

with initial data 5° = r\ — r 2 . Writing a similar energy estimate to the 
one above yields formally 

|(l|A^||i 2 + ||V^||i 2 )+^||V ft A^||i 2 

= - / V^S ■ V h (A hri )A h S dx. 

Jm. 2 

Then we simply write, by Holder's inequality followed by Gagliardo- 
Nirenberg's inequality 



/ Vi5-V h (A h n)A h 5dx <\\ViS\\ L 4V h A h n\\ L 2\\A h S\\ L 4 

<C||V^||| 2 ||A^||J 2 ||V /l A ft r 1 |U2||A^||J 2 ||V ft A ft 5||| 2 . 
This implies that 

/ V£5-V h (A h n)A h 5dx <=\\V h A h 5\\ 2 L2 + \\V h 5\\ 2 L2 
Jm 2 Q 

a rll2 ||V7 A ». 112 



+C A /J||A h< y||i 9 ||V fc A h r 1 ||i a . 



Finally Gronwall's inequality allows to obtain 

l|A^(t)||| 2 + \\V h 5(t)\\l 2 + ±f \\V h A h 8{t')f L . dt' 

Q Jo 

|A^°||| 2 + ||V^°||i 2 ) exp (c^J \\V h A hri (t')f L2 dt' + Ct). 



This allows to conclude to stability, hence uniqueness for the limit 
system (leaving the usual regularization procedure to make the above 
arguments rigorous to the reader) provided the initial datum enjoys 
the extra regularity stated in Theorem 12.11 
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